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As  a  result  of  increasing  frequency  and  intensity  of  tropical  cyclones,  an  accurate  forecasting  of 
cyclone  evolution  and  ocean  response  is  becoming  even  more  important  to  reduce  threats  to  lives 
and  property  in  coastal  regions.  To  improve  predictions,  accurate  evaluation  of  the  air-sea 
momentum  exchange  is  required.  Using  current  observations  recorded  during  a  major  tropical 
cyclone,  we  have  estimated  this  momentum  transfer  from  the  ocean  side  of  the  air-sea 
interface,  and  we  discuss  it  in  terms  of  the  drag  coefficient.  For  winds  between  20  and  48  meters 
per  second,  this  coefficient  initially  increases  and  peaks  at  winds  of  about  32  meters  per 
second  before  decreasing. 


The  air-sea  momentum  exchange  under  a 
tropical  cyclone  determines  the  oceanic 
response  to  its  winds.  An  accurate  estima¬ 
tion  of  the  exchange  hence  is  required  for 
correctly  forecasting  storm  track  and  intensi¬ 
ty,  as  well  as  for  accurately  predicting  storm 
surges,  ocean  currents,  and  waves,  and  for 
making  hurricane  risk  assessments,  particularly 
because  the  frequency  and  intensity  of  tropical 
cyclones  are  reportedly  increasing  (/).  Con¬ 
ventional  methods  of  determining  the  air-sea 
momentum  transfer  are  from  the  atmospheric 
side  of  the  interface  and  are  based  on  mea¬ 
surements  of  wind  profiles  near  the  ocean  sur¬ 
face  or  of  wind  turbulence  (2).  These  methods, 
however,  cannot  work  reliably  for  the  extreme 
high-wind  conditions  of  a  major  tropical  cy¬ 
clone,  because  wind  measurements  near  the 
ocean  interface  are  inevitably  contaminated 
from  intense  breaking-wave  activities  and  sea 
spray  (J). 

Most  available  field  measurements  of  sur¬ 
face  stress  are  for  wind  speeds  ranging  from  4  to 
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20jns  '.  The  techniques  employed  estimate  the 
interface  momentum  exchange  as  a  wind  stress,  rs, 
given  by 

Tv  =  P^CojW]-  ( I ) 

where  pajr  is  the  air  density,  CD  is  the  drag 
coefficient,  and  |W|  is  the  velocity  magnitude  of 
the  wind  at  a  reference  height  (usually  10  m). 

In  open  ocean  conditions,  as  the  wind 
becomes  stronger  the  range  of  surface  waves, 
which  travel  at  slower  speeds  than  the  wind, 
widens.  Accordingly,  CD  increases  as  wind  speed 
increases  (J).  This  pattern  of  variability  in  CD 
is  commonly  used  in  wave,  surge,  and  circula¬ 
tion  numerical  models,  which  incorporate  a 
monotonic  increase  of  the  drag  coefficient  with 
wind  speed  even  under  extreme  wind  conditions. 
Such  an  assumption  may  result  in  an  over¬ 
estimation  of  the  air-sea  momentum  transfer  and 
lead  to  unrealistic  model  predictions  or.  perhaps 
even  worse,  model  tuning  using  bad  physical 
justifications. 

Recently,  limited  meteorological  observa¬ 
tions,  theory,  and  experiments  (3-8)  suggest  a 
reduction  in  the  drag  coefficient  at  high  wind 
speeds,  generally  greater  than  30  ms'1,  All  these 
efforts  to  evaluate  the  air-sea  momentum  ex¬ 
change  have  been  made  through  meteorological 
measurements,  or  theories  and  models  developed 


for  the  atmospheric  boundaiy  layer  (“top-down” 
determination).  Estimation  of  this  exchange 
can  also  be  done  from  the  other  side  of  the 
interface  with  full  water-column  ocean  current 
observations.  This  method  directly  determines 
the  stress  at  the  interface  based  on  ocean 
currents  and  provides  a  very  reliable  and 
accurate  direct  determination  of  the  air-sea 
momentum  transfer  under  high  winds  (3). 
However,  this  type  of  “bottom-up”  determina¬ 
tion  imposes  the  almost  impossible  require¬ 
ments  of  deploying  instruments  in  the  ocean 
directly  under  the  highly  unpredictable  path  of 
a  major  tropical  cyclone  and  then  having  the 
instruments  survive  the  enormous  forces  gen¬ 
erated  by  the  cyclone. 

On  15  September  2004,  the  center  of 
Hurricane  Ivan  (at  category  4  strength)  passed 
directly  over  six  current  and  wave/tide  gauge 
moorings  on  the  outer  continental  shelf  in  the 
northeastern  Gulf  of  Mexico  (9-12)  (Fig.  1). 
The  instrumentation  survived  and  measured 
nearly  full  ocean  current  velocity  profiles 
directly  under  the  cyclone.  The  resulting  data 
set  allows  direct  “bottom-up”  determination  of 
the  air-sea  momentum  exchange  under  extreme 
wind  conditions. 

Extreme  winds  of  a  tropical  cyclone  elicit  a 
four-stage  response  when  passing  over  ocean 
waters  (9, 13).  The  first  three  stages  make  up  the 
"forced  stage”  response,  whereas  the  fourth  stage 
is  the  “relaxation  stage."  Over  the  open  ocean, 
the  forced  stage  response  is  primarily  baroclinic 
(depth-dependent),  with  a  weak  barotropic 
(depth-independent)  response  consisting  of  a 
trough  in  sea  surface  height  and  an  associated 
geostrophic  'current  that  are  set  up  almost 
instantly  (13).  Over  the  continental  shelf,  how¬ 
ever,  the  forced  response  is  quite  different  (9).  It 
consists  of  a  strong  barotropic  component  that 
is  not  geostrophically  balanced  and  a  much 
weaker  baroclinic  response.  Thus,  the  directly 
forced  response  on  the  continental  shelf  should 
be  well  described  to  the  first  order  by  the  linear 
time-dependent  depth-integrated  horizontal  mo¬ 
mentum  equations.  Scaling  analysis  of  current 
velocity  measurements  recorded  in  the  northeast¬ 
ern  Gulf  of  Mexico  (Fig.  1 )  under  Humcane  Ivan 
supports  this  assumption,  especially  for  the  along- 
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Fig.  1.  Regional  map 
and  instrument  locations. 
Bathymetry  (in  meters) 
is  shown  for  the  north¬ 
eastern  Gulf  of  Mexico. 
Acoustic  Doppler  current 
profiler  (ADCP)  moorings 
are  denoted  by  black  dots 
(Ml  to  M6).  The  yellow 
triangle  is  Buoy  42040, 
National  Data  Buoy  Cen¬ 
ter.  Hurricane  Ivan's  path 
is  indicated  by  the  dashed 
red  line.  Two  inserts  in  the 
expanded  map  are  a 
colorized  infrared  image 
of  Hurricane  Ivan  taken 
by  the  National  Oceanic 
and  Atmospheric  Admin¬ 
istration's  Geostationary 
Operational  Environmen¬ 
tal  Satellite-12  satellite 
(lower  right)  and  a  Barny 
mooring  containing  the 
ADCP  current  meter  (up¬ 
per  right). 


Fig.  2.  Drag  coefficient 
for  several  resistant  coef¬ 
ficients.  C0  is  shown  as  a 
function  of  the  wind  speed 
at  10  m  (W10)  for  several 
different  resistance  coeffi¬ 
cients  r  (open  black  circles 
—  r  —  0.001  cm  s'1;  green 
triangles  -  r  =  0.0505  cm 
s'1;  blue  asterisks  -  r  = 
0.1  cm  s'1).  The  solid  lines 
represent  quadratic  curves 
fitted  to  the  evaluated  C0 
separately  for  each  r. 


shelf  momentum  [more  discussion  concerning 
the  along-shelf  momentum  balance  is  in  (/4)] 
given  by 


dU_ 

dt 


rU 

H 


(2) 


where  p  is  a  reference  density  ( 1 025  kg  m'3),/ 
is  the  Coriolis  parameter  (0.71  x  10'4  s*1),  U 
and  Fare  the  depth-integrated  along-shelf  and 
cross-shelf  velocity  components,  H  is  the  water 
depth,  r  is  a  constant  resistance  coefficient  at 


the  sea  floor,  and  tsx  is  the  along-shelf  wind 
stress. 

The  simplified  dynamics  given  by  Eq.  2  may 
not  hold  under  all  conditions.  For  instance,  in  this 
case,  the  momentum  balance  breaks  down  once 
the  eye  of  the  storm  has  passed  over  the  moorings 
and  the  wind  vectors  rapidly  rotate  180°.  When 
this  happens,  the  water  column  cannot  respond 
and  change  direction  as  rapidly  as  the  wind,  and 
these  dynamics  no  longer  hold.  Additionally,  the 
momentum  balance  may  not  hold  for  a  rapidly 
moving  tropical  cyclone  or  for  a  cyclone  with  an 


erratic  path  for  similar  reasons.  However,  many 
cyclones  move  slowly  and  steadily  enough 
over  the  shelf  for  such  an  approximation  to 
hold. 

The  air-sea  momentum  transfer,  expressed  in 
terms  of  tw,  is  estimated  from  Eq.  2  using  our 
ocean  current  observations.  In  atmospheric 
studies,  this  exchange  is  commonly  discussed  in 
terms  of  the  drag  coefficient  (CD).  Wind  stress 
(Tix)  can  be  also  defined  as  t*,  =  pmCD\W\lVx, 
where  Wx  is  the  along-shelf  wind  velocity 
component.  Thus,  for  comparison  with  other 
studies,  we  can  determine  CD  by  inserting  the 
formula  for  r„  into  Eq.  2. 

r  -  pH  (eU  nr  rU\ 

D  Pjwin^a  f  +  77 J  (3) 

Results  produced  from  evaluation  of  Eq.  3 
are  shown  in  Figs.  2  and  3.  The  drag  coeffi¬ 
cient  was  derived  for  several  different  values 
of  the  resistance  coefficient  (r)  ranging  between 
0.001  cm  s'1  and  0.1  cm  s'1.  These  values  of  r 
are  similar  to  those  used  in  other  studies  inves¬ 
tigating  shelf  currents  (15-18).  Figure  2  displays 
estimates  of  CD  for  minimum  (0.001  cm  s"1), 
mean  (0.0505  cm  s'1),  and  maximum  (0. 1  cm  s'1) 
values  of  r  used  in  Eq.  3,  whereas  Fig.  3  shows 
results  for  r  —  0.02  cm  s  1  .that  is,  the  resistance 
coefficient  that  was  determined  from  observa¬ 
tions  for  the  northeastern  shelf  of  the  Gulf  of 
Mexico  (15).  The  results,  especially  for  the  wind 
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Fig.  3.  Drag  coefficient  as 
a  function  of  wind  speed.  CD 
is  shown  for  an  observation- 
based  resistance  coefficient, 
r  =  0.02  cm  s_1.  The  red 
open  circles  are  the  eval¬ 
uated  C0  from  the  current 
and  wind  observations,  the 
solid  red  line  is  a  fitted 
quadratic  curve  to  the  C0 
estimates,  and  the  red 
dashed  lines  are  the  95% 
confidence  limits  for  this 
quadratic  curve.  The  black 
dotted  tines  represent  the 
window  for  C0  reported  in 
(6),  whereas  the  blue  dots 
represent  CD  reported  in  (4). 


Wind  spend  {W10,  m  *'') 


speeds  below  30  ms"1,®  somewhat  noisy  as  a 
result  of  measurement  uncertainty  and  the  need 
to  calculate  a  velocity  derivative,  which  tends  to 
enhance  noise.  However,  they  consistently 
show  a  decreasing  trend  of  Cp  for  wind  speeds 
greater  than  32  m  s_1,  the  lower  threshold  for  a 
category  1  hurricane  on  the  Saffir-Simpson  Scale. 
It  is  also  apparent  that  the  CD  values  are  weakly 
dependent  on  the  choice  of  the  resistance  co¬ 
efficient  and  are  larger  for  increasing  values 
of  r.  The  drag  coefficient  estimates  evaluated 
for  r  =  0.1  cms'1  are,  on  average,  20%  greater 
than  those  calculated  for  r  =  0.001  cm  s-1  from 
Eq.  3. 

To  produce  the  best  representation  of  Cp  for 
each  r,  a  second-order  curve  (a  function  of  the 
wind  speed)  was  fitted  by  a  least-squares 
technique  to  all  estimated  values  of  Cp.  The 
curves  are  displayed  in  Figs.  2  and  3.  Addition¬ 
ally,  the  95%  confidence  limits  for  the  fitted 
curve  are  shown  in  Fig.  3.  The  pattern  of  the 
relationship  between  Cp  and  the  wind  speed  is 
robust,  but  the  curve  coefficients  are  determined 
by  the  value  chosen  for  r  in  Eq.  3.  However,  ail 
curves  clearly  show  an  initial  increase  of  the  drag 
coefficient  and  monotonic  decrease  as  found  by 
recent  studies  (3-8)  after  reaching  a  maximum 
value  at  ~32  m  s-1.  Some  of  these  studies  (3,  19) 
imply  that  the  decreasing  drag  at  high  winds 
seems  to  be  related  to  the  spray,  foam,  and 
bubbles  from  breaking  waves  that  reduce  the 
drag  and  allow  the  hurricane  to  slip  over  the  sea. 

With  the  nearly  full  water-column  ocean  cur¬ 
rent  measurements,  the  only  unknown  term  left 
in  the  simplified  equation  of  motion  is  the  wind 
stress.  Thus,  the  behavior  of  the  drag  coefficient 
(Cp)  can  easily  be  estimated  for  a  range  of  strong 
winds.  Despite  the  fact  that  the  drag  coefficient  is 
evaluated  differently  here,  estimates  of  Cp 
determined  “bottom-up”  reasonably  replicate 
the  values  determined  “top-down”  in  recent 
studies  (3-7),  Results  from  our  research  show 
that  Cp  peaks  at  a  wind  speed  near  32  m  s-1  and- 


then  steadily  decreases  as  the  wind  speed 
continues  to  rise.  Our  values  for  CD  are  in  a 
range  of  Cp  values  found  using  meteorological 
observations  (4)  for  wind  speeds  greater  than 
32  m  s'1  but  are  higher  for  lower  wind  speeds. 
These  differences  may  be  attributed  to  uncertain¬ 
ties  in  the  wind  measurements  and  the  applica¬ 
bility  of  the  simplified  ocean  dynamics  at  the 
lower  wind  speeds. 
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